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Abstract. The angular distribution of the emitted electrons, following the two- 
photon ionization of the hydrogen- like ions, is studied within the framework of second 
order perturbation theory and the Dirac equation. Using a density matrix approach, we 
have investigated the effects which arise from the polarization of the incoming light as 
well as from the higher multipoles in the expansion of the electron-photon interaction. 
For medium- and high-Z ions, in particular, the non-dipole contributions give rise to 
a significant change in the angular distribution of the emitted electrons, if compared 
with the electric-dipole approximation. This includes a strong forward emission while, 
in dipole approxmation, the electron emission always occurs symmetric with respect 
to the plane which is perpendicular to the photon beam. Detailed computations for 
the dependence of the photoelectron angular distributions on the polarization of the 
incident light are carried out for the ionization of H, Xe^^+, and U^^+ (hydrogen- like) 
ions. 



1. Introduction 

During the last decades, the multi-photon ionization of atoms and ions has been 
widely studied, both experimentally and theoretically. While, however, the majority 
of experiments were first of all concerned with the multi-photon ionization of complex 
atoms, most theoretical investigations instead dealt with the ionization (and excitation) 
of the much simpler hydrogen-like and helium-like systems. For atomic hydrogen, in 
contrast, multi-photon experiments have been carried out only recently (Wolff et al 
1988, Rottke et al 1990, Antoine et al 1996) because of the former lack of sufficiently 
intensive (and coherent) light sources in the UV and EUV region. With the recent 
progress in the setup of intensive light sources in the EUV and x-ray domain, such as the 
fourth-generation synchrotron facilities or variously proposed free-electron lasers, two- 
and multi-photon studies on the ionization of inner-shell electrons are now becoming 
more likely to be carried out in the future, including case studies on medium-Z 
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and high-Z hydrogen-like ions (Kornberg et al 2002). With increasing charge (and 
intensity of the hght), of course, relativistic effects will become important and have 
been investigated in the past for the two-photon excitation and decay (Goldman and 
Drake 1981, Szymanowski et al 1997, Santos et al 2001) as well as ionization (Koval et 
al 2003) of hydrogen-like ions. So far, however, all of these studies were focused on the 
total (excitation or decay) rates and ionization cross sections while, to the best of our 
knowledge, no attempts have been made to analyze the effects of relativity on angular 
resolved studies. 

In this contribution, we explore the angular distribution of the electrons following 
the two-photon ionization of hydrogen-like ions. Second-order perturbation theory, 
based on Dirac's equation, is applied to calculate the two-photon amplitudes including 
the full (relativistic) electron-photon interaction. The angular distribution of the photo- 
electrons are then derived by means of the density matrix theory which has been found 
appropriate for most collision and ionization processes and, in particular, for angular- 
dependent studies (Laplanche et al 1986). Since, however, the basic concepts of the 
density matrix theory has been presented elsewhere at various places (Blum 1981, 
Balashov et al 2000), we will restrict ourselves to rather a short account of this theory 
in Subsection 12. 1[ Apart from a few basic relations, here we only show how the angular 
distribution of the electrons can be traced back to the two-photon transition amplitudes. 
The evaluation of these amplitudes in second-order perturbation theory and by means of 
Coulomb-Green's functions are discussed later in Subsections 12. 2l and l2.3[ and including 
the full decomposition of the photon field in terms of its multipole components in 
Subsection 12. 4[ Using such a decomposition, we have calculated the electron angular 
distributions for the two-photon ionization of the Is ground state of hydrogen (H) as 
well as hydrogen-like xenon (Xe^^"*") and uranium (U^^^). By comparing the angular 
distributions for different nuclear charges Z, we were able to analyze both, the effects 
of the polarization of the — incoming — light and the contributions from higher (i.e. 
non-dipole) multipoles in the decomposition of the electron-photon interaction. These 
results are displayed in Section [3] and clearly show that, with increasing charge Z, the 
higher multipole components lead to a strong emission in forward direction (i.e. parallel 
to the propagation of the light), while the electric-dipole approximation alone gives rise 
to a symmetric electron emission around the polar angle 9 90°, similar as obtained 
by nonrelativistic computations (Zernik 1964, Lambropoulus 1972, Arnous et al 1973). 
Finally, a brief summary on the two-photon ionization of medium and high-Z ions is 
given in Section |H 

2. Theory 

2.1. Density matrix approach 

Within the density matrix theory, the state of a physical system is described in terms of 
so-called statistical (or density) operators (Fano 1957). These operators can be consi- 
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dered to represent, for instance, an ensemble of systems which are — altogether — 
in either a pure quantum state or in a mixture of different states with any degree of 
coherence. Then, the basic idea of the density matrix formalism is to accompany such 
an ensemble through the collision process, starting from a well defined 'initial' state and 
by passing through one or, possibly, several intermediate states until the 'final' state of 
the collision process is attained. 

In the two-photon ionization of hydrogen-like ions, the 'initial' state of the (com- 
bined) system 'ion plus photons' is given by the bound electron \nbjbfJ'b) o-nd the two 
incoming photons, if we assume a zero nuclear spin J = 0. For the sake of simplicity, 
we also restrict our treatment to the case that both photons will have equal momentum: 
ki = k2 = k, while the spin states of the photons may still differ from each other and 
are characterized in terms of the helicity parameters Ai, A2 = ± 1 (i.e. by means of 
their spin projections onto the direction of propagation k). Of course, the case of equal 
photon momenta k correspond to the most frequent experimental setup of the two- 
photon ionization of atoms and ions using, for instance, lasers or synchrotron radiation 
sources. With these assumptions in mind, the initial spin state of the overall system 
is determined by the direct product of the statistical operators of the ion and the two 
incident photons 

pi= pb®P-i® Pi (1) 

or, explicitly, in a representation of the density matrix in terms of the individual 
momenta by 

{ubjbPb, kAi, kAs I Pi I Ubjbp'b, kA'i, kX'^) 

= {nbjbPb\pb\nbjbPb) (kAilp^lkA'i) (kA2|p^|kA2) . (2) 

In the 'final' state of the ionization, after the electron has left the nucleus, we just 
have a free electron with asymptotic momentum p and spin projection (as well as 
the bare residual ion with nuclear charge Z). Therefore, the final spin state is described 
by the statistical operator of the emitted (free) electron pe which, in the framework of 
the density matrix theory, can be obtained from the initial-state density operator pi 
owing to the relation (Blum 1981, Balashov 2001) 

Pf = Pe = RpiR^ ■ (3) 

In this simple relation, R is called the transition operator and must describe the 
interaction of the electron with the (two photons of the) radiation field. Of course, 
the particular form of the transition operator R depends on the framework in which we 
describe the coupling of the radiation field to the atom. As appropriate for high-Z ions, 
below we will always refer to a relativistic treatment of the electron-photon interaction, 
based on Dirac's equation and the minimal coupling of the radiation field (Berestetskii 
et al 1971). 

Instead of applying Eq. ([3]), in practise, it is often more convenient to rewrite the 
statistical operators in a matrix representation. Using, for example, the initial spin 
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Figure 1. Geometry of the two-photon ionization. The photoelectron is emitted 
along the unit vector p = {6, (p) where 6 is the (polar) angle between the incident 
photon momenta k (chosen as the z-axis) and the electron momentum p. Moreover, 
the (azimuthal) angle (f) defines the angle of p with respect to the x-z plane which, in 
the case of linearly-polarized light, contains the polarization vector u. 



density matrix (l2l), we easily obtain the density matrix of the (finally) emitted electron 
by 

(pm J pe I P"^s) = ^ X] (nfojbpb I pb I rifojb/i'f,) (kAi I I kA'i) (kA2 | p^ | kAg) 

X Mf,p(m„p6, Ai, A2) M6*p(m;,p'b, A'l, A2), (4) 
where used is made of the abbreviation 

Mbp(m^,Pfe, Ai, A2) = (pm^ R kAi, kA2, rifejfePfe) (5) 



in order to represent the transition amplitudes for the two-photon ionization. The final- 
state density matrix (jl]) still contains the complete information about the ionization 
process (i.e. the properties of the bare ion and the electron) and, thus, can be used 
to derive all the observable properties of the photoelectrons. Obviously, however, the 
outcome of some considered experiment will depend on the particular setup and the 
capability of the detectors for resolving the individual properties of the particles. In the 
density matrix theory, this setup of the experiment is typically described in terms of a 
(so-called) detector operator P which characterizes the detector system as a whole. In 
fact, this detector operator can be considered to project out all those quantum states of 
the final-state system which leads to a 'count' at the detectors; in the language of the 
density matrix, therefore, the probability for an 'event' at the detector is simply given 
by the trace of the detector operator with the density matrix: W = Tr(P p) . 

To determine, for instance, the angular distribution of the emitted (photo-) elec- 
trons, we may assume a detector operator in a given direction p = {9, 0) [cf. Figure [1] 
which is insensitive to the polarization of the electrons 

P = ^ Ipm^) (pm^l , (6) 
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i.e. a projection operator along p and including a summation over the spin state 
of the electrons. From this operator, and by taking the trace over the product (P p) 
with the final-state density matrix (j4j), we obtain immediately the electron angular 
distribution in the form 

W{p) = TriPpf) = E (kAilp.lkA;) (kA^ kA'^) 

Mtms A1AJA2A2 

X Mt,p{ms, Hb, Ai, A2) M^*p(ms, /if,, A'^, A'a) (7) 

where, for the sake of simplicity, we have assumed that the hydrogen-like ion is initially 
unpolarized. Apart from this additional assumption, however, Eq. ((Tj) still represents 
the general form of the electron angular distribution for the process of the two-photon 
ionization of hydrogen-like ions. As seen from this equation, the emission of the 
photoelectron will depend on the spin state of the incident photons, defined by the 
photon density matrices (kA|p^|kA'). For any further evaluation of this distribution 
function, therefore, we shall first specify these density matrices or, in other words, the 
polarization of the incoming light. For example, if both photons are unpolarized, the 
(two) photon density matrices simply reduce to a constant 1/2, (kA | \ kA') = 5aa'/2 
[cf. Appendix, Eq. IA.2j and leads us to the well-known angular distribution 

^'^"^(P) = ^7^-^ E E |M,pK,/.,,Ai,A2)r- (8) 

^ ^ fJ-b-rris A1A2 

For many (modern) light sources such as lasers or synchrotron radiation, it is not 
very practical to consider only unpolarized light from the very beginning. In general, 
instead, the angular distribution of the emitted electrons will depend both on the type 
as well as the degree of the polarization of the incident light. For circularly polarized 
light with degree Pc, for instance, the photon density matrix from Eq. ([7]) becomes 
(kA I p^ I kA') = (1 + XPc) 5aa'/2 and, hence, give rise to the angular distribution 

W^^'iP) = E (1 + ^1^^) (1 + ^2^^) |M,p(m„/x,,Ai,A2)|' , (9) 

while, for linearly polarized light along the and with a polarization degree Pl, 

the photon density matrix is (kA | py \ kA') = 5aa'/2 + (1 "~ 5\\')PlI'^- If we evaluate 
Eq. (171) again with this latter density matrix, we obtain the angular distribution 



^ ^ fitms \ A1A2 



Mbp{ms,f^b, Ai, A2) 1^ 



Pi 



{■ms,fib, Ai, A2) 

Ai A2 



2N 

^ (10) 



for the electrons as emitted in the two-photon ionization of hydrogen-like ions with 
linearly polarized light. 
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2.2. Two-photon transition amplitude in second- order perturbation theory 

For any further analysis of the electron angular distributions, following the two- 
photon ionization of a hydrogen-like ion, we need to calculate the transition amplitude 
Mbp{ms, ^b, Xi, \2) as seen from Eqs. (!8!)- (fT0l) . This amplitude describes a bound-free 
transition of the electron under the (simultaneous) absorption of two photons. For 
a moderate intensity of the photon field, of course, this amplitude is most simply 
calculated by means of second-order perturbation theory (Laplanche et al 1976) 

M (m n \ \\ Sf ^'^P'"- I" ""^^ e^''^! \a UA, e''''\ipn,j,t.,) , . 

Mbp[ms, Hb, ■^2) - 2L u' B B ' ^-'--'-^ 

where the transition operator a e***^^ describes the (relativistic) electron-photon inter- 
action, the unit vector ua the polarization of the photons, and where the summation 
runs over the complete one-particle spectrum. From the energy conservation, moreover, 
it follows immediately that the energies of the initial bound state, Eb, and the final 
continuum state, Ef, are related to each other by Ef = Eb + 2£'^, owing to the 
energy of the incoming photons, E^ = hck. Although known for a long time, the 
relativistic form of the transition amplitude ffTTl) has been used only recently in studying 
multi-photon ionization processes and, in particular, in order to calculate the total 
ionization cross sections along the hydrogen isoelectronic sequence (Koval et al 2003). 
In such a relativistic description of the transition amplitude f fTTj) . the initial state 
'^Hiii/ii, (r) = (r I ribjbfJ^b) and the final state ippmA''^) = (''I P^s) are the (analytically) 
well-known solutions of the Dirac Hamiltonian for a bound and continuum electron, 
respectively (Berestetskii et al 1971). 

As seen from Eq. ( [TT]) . the evaluation of the transition amplitude requires a 
summation over the discrete (bound) states as well as an integration over the continuum 
of the Dirac Hamiltonian, {tp,y,E,y). In fact, such a 'summation' over the complete 
spectrum is difficult to carry out explicitly since, in particular the integration over 
the continuum requires the calculation of free-free transitions. In the past, therefore, 
this summation has often been restricted to some small — discrete — basis, assuming 
that the contribution from the continuum is negligible. In practise, however, such a 
limitation seems justified only to estimate the behaviour of the cross sections near the 
resonances where the ion is rather likely excited by the first photon into some — real — 
intermeditate state of the ion from which it is later ionized by means of a second photon. 
In the non-resonant region of the photon energies, in contrast, the integration over the 
continuum may give rise to a rather remarkable contribution to the total cross section 
and, hence, has to be carried out. Apart from a direct summation over the continuum 
states, however, it is often more favorable to apply Green's functions, at least if these 
functions can be generated efficiently. For hydrogen-like ions, for example, such Green's 
functions are known analytically, both in the nonrelativistic as well as relativistic theory 
(Swainson and Drake 1991). 
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2.3. Green's function approach 

As usual, Green's functions are defined as solutions to some inhomogeneous (differential) 
equation 

'E-H]GE{r,v') = 5{v-v') (12) 



where, in our present investigation, H refers the Dirac Hamiltonian and E denotes the 
energy of the atom or ion. For realistic systems, of course, such Green's functions are 
not easy to obtain, even if only approximate solutions are needed. However, a formal 
solution is given by (Morse and Feshbach 1953) 

G,(...-)^^ i^-';»<yi . (13, 

including a summation (integration) over the complete spectrum (of K) as discussed in 
the previous section. In the two-photon transition amplitude (11 II) . therefore, we may 
simply replace this summation by the corresponding Green's function 

Mfcp(m^,/Xb, Ai, Aa) = (^V^pm,(r) «UAie*''''G'ij^+ij^(r, r') auAje*''''' V„^K^^Jr')^ . (14) 

For hydrogen-like ions, the Coulomb-Green's functions from Eq. f|T2|) are known 
analytically today in terms of (various) special functions from mathematical physics 
and, in particular, in terms of the confluent hypergeometric function iFi(a, 6, z). Here, 
we will not display these functions explicitly but refer the reader instead to the literature 
(Swainson and Drake 1991, Koval and Fritzsche 2003). For the further evaluation of 
the transition amplitudes f[T^ let us note only that, also for the one-particle Dirac 
Hamiltonian, the Coulomb-Green's function can be decomposed into a radial and an 
angular part 

. \ sr^ ( gk^ir,r')n^^if)niUf') -igk'Jr,r') Q^^r) ^l.^ir') \ , , 
GE[r,r) = — y [ ,(15) 

^ \ ig|^(r,r') 0_.™(r) l^L(r') g|f(r,r') 0_.„(f) ^l,^(f') ) 

where the ^f^^if) denote standard Dirac spinors and where the radial Green's function 
is given in terms of four components g^^'(r, r') with T = L,S refering to the large and 
small components of the associated (relativistic) wave functions. The computation of 
the radial Green's function for hydrogen-like ions has been described and implemented 
previously into the Greens library (Koval and Fritzsche 2003); this code has been 
used also for the computation of all transition amplitudes and (angle-differential) cross 
sections as shown and discussed below. 



2.4- Exact relativistic formulation of the two-photon amplitude 

Eq. (1141) displays the two-photon transition amplitude in terms of the (relativistic) wave 
and Green's functions of hydrogen-like ions. For the further evaluation of this amplitude, 
we need to decompose both, the photon as well as the free-electron wave functions 
into partial waves in order to make later use of the techniques of Racah's algebra. 
As discussed previously for the capture of electrons by bare, high-Z ions (Surzhykov 
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et al 2002), we first have to decide about a proper quantization axis (z-axis) for tfiis 
decomposition, depending — of course — on tfie particular process under consideration. 
For tlie pliotoionization of atoms, the only really prefared direction of the overall system 
is given by the photon momenta ki = k2 = k which we adopt as the quantization axis 
below. Then, the multipole expansion of the radiation field reads as 

oo 

u.e^'^r = u,e"=^ = V2^Y.'' t^]'^' ( ^lx^ + ^AAg ) , (16) 

L=l 

where [L] = (2L + 1) and the standard notation A^^^™"* is used for the electric and 
magnetic multipole fields, respectively. Each of these multipoles can be expressed in 
terms of the spherical Bessel functions jhikr) and the vector spherical harmonics T|^^ 
of rank L as (Rose 1957): 

AS=J.(A:r)T^„ 



Ait = J.-i(fcr) ^^^TiVi-JL+i(A:r) ^^Ti^L+i- (17) 

Using the expressions f|T6|) and f|T7j) for the photon field, we can rewrite the two- 
photon transition amplitude ( |T4l) in terms of its electric-magnetic components 



oo 

L,L'=1 AA' 

X (V-pm. \ajA{kr) Tl]j, GE,+E,{r,r') a jA'{kr') | ^n,.,^,> , (18) 

where the coefficients are defined as 

1 if A = L 



e^A= ^Virr ifA = ^-i . (19) 



-^^Vra ifA = L + l 
As seen from the expansion f|T8l) . we can distinguish between different multipole compo- 
nents such as eIeI, eImI, e1e2, and others owing to the symmetries of the two vector 
spherical harmonics, i.e. due to the particular combination of the summation indices 
L, L', A, A' in this expansion. In the second line of (fT8l) . however, the — electro-magnetic 
— multipole matrix elements still contain the wave function tppms (r) of the free electron 
with well-defined asymptotic momentum p. In another expansion, therefore, we have to 
decompose it into partial waves to allow for a further simplification of the two-photon 
transition amplitude (ITS!) . Again, also the expansion of the free-electron wave will 
depend on the choice of the quantization axis and requires — by using a quantization 
along the photon momentum — that we have to carry out a rotation of the space part 
of the electron wavefunction from the z-direction into the p-direction 

V'pm,(r) = 47r ^ i'/ e'^^^f {ifj^j - rris I /2ms \ jfjJ'f) 

I igL,Mf!-.„ (ft) ^ 
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where the summation runs over all partial waves kj = ±1,±2..., i.e. over all possible 
values of the Dirac angular momentum quantum number Kf = ±(jj + 1/2) for 
I J = jf ± 1/2. In this notation, the (nonrelativistic angular) momentum represents 
the parity of the partial waves and A^^ is the Coulomb phase shift. Moreover, as seen 
from expression (l20l) . the partial waves 

separate into a radial and an angular parts, where the two radial functions 

g|,(r) = PE.{r), gL(r) = g^«(r) 

are often called the large and small components and the corresponding angular parts 
^Kf,,f{n) = \lfjff^f) = J2m,ms {^f ™i 1/2 I jf /i/> Yifn^S^) Xi/2 and n.^f^^in) = 
\lfjflJ^f) = Y.mims {^f 1/2 I jf At/) ^«s„^(n) Xi/2 m, are the standard Dirac spin- 
angular functions. 

Using the partial-wave decomposition fl2T|) for the free-electron wave function and 
a similar expansion f|T5|) for the Green's functions, we now can carry out the angular 
integration in the transition amplitude ( |T8l) analytically 

M,K,/.,,Ai,A2) = 87r2 E E E P^' e'^^f 

LkL'k' Kffif KmTT' 

Kflj fif \a T kI m) (kI m | a T^?^, | k^/^ /if, 

X U 1% (/€/, K, Kb) YiL (n) (22) 

where, apart from the Clebsch-Gordan coefficient {ijfJ^f — mg l/2ms | jffJ^f) and some 
constant factors, the angular part of the amplitude is given in terms of the matrix 
elements of the rank L spherical tensor cxT^^ = [Ya ■ These matrix elements can 

be simplified to (Balashov et al 2000) 



jbfJ'b) 



If 1/2 Jb 
'4^0, AO 1/^0) {A 1 L } , (23) 

e 1/2 Ja 

by using a proper decomposition in terms of the orbital and spin subspaces. The radial 
part of the transition amplitude (ITSjl is contained in (l22il in the (two-dimensional) 
integrals 

[/JJ;(k^, «:,«:,) = J gi^,^(r) jA(A;r) gif;^^Jr,r') jA,(A;r') gj,^(r') rfr rfr', (24) 

which combines the various (large and small) components of the bound state, the Green's 
function as well as from the free-electron wave. In this notation, again, T = L, S and 
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a superscript T refers to the conjugate of T, i.e. T = L for T = S and vice versa. 
In contrast to the angular integrals (|23|) . the radial integrals have to be computed 
numerically. In the present work, all the required integrals for the two-photon transition 
amplitudes (122|) are calculated by using the Greens (Koval and Fritzsche, 2003) and 
Racah (Fritzsche et al 2001) programs. 

2.5. Electric dipole approximation 

The transition amplitude fl22|) still describes the full interaction between the electron 
and photon fields. With the explicit summation over all the multipoles of the photon 
field (fT6|) . it includes the so-called retardation effects or non-dipole contributions. In 
practise, however, the contributions from the higher multipoles decreases very rapidly 
with L and may therefore be neglected; in fact, the computation of these contributions 
also become rather tedious because of difficulties with a stable procedure for the two- 
dimensional radial integrals ( l24l) . In many cases, therefore, it seems justified to restrict 
the summation in (!22l) to just the (dominant) electric dipole term with L = 1 and 
A = L± 1. This 'dipole approximation' is valid if the photon wave length is much larger 
than the size of the atom, i.e. ka^ -C 1 where Oq is the Bohr radius. For the two-photon 
ionization, this condition is well satisfied for most light ions with, say, Z < 30 and for 
photon energies below of the one-photon ionization threshold. 

From the general form (1221) of the ionization amplitude, the electric dipole appro- 
ximation is obtained by taking L = L' = 1 and A, A' = 0, 2 which — owing to the dipole 
selection rules — then also restricts the summation over k/, i.e. the allowed partial waves 
for the free electron. For the ii"-shell ionization with (completely) circularly polarized 
light, for instance, the final-state electron can only escape in the ^3/2 or ^5/2 states. 
And, as seen from Eq. (l22l) . the dipole transition amplitude is then indeed defined by 
the (second-rank) spherical harmonic, Mfep(ms, /i^. A, A) oc F2, /ij,-ms+2A(p) ~ sin^iO) 
which (together with eq. [9]) leads us to the well-known angular distribution 



of the photoelectrons (Lambropoulos 1972, Arnous et al 1973). As expected from the 
axial symmetry of the overall system 'ion plus photons', the angular distribution (!25ll 
only depends on 6 but not on the azimuthal angle (f). For linearly polarized light, in 
contrast, a reaction plane is naturally defined by the photon momentum k and the pola- 
rization vector u and, hence, the axial symmetry is broken. For a linear polarization 
of the incident light, therefore, the angular distribution will depend on both, the polar 
and azimuthal angle and is given by (Zernik 1964, Lambropoulos 1972) 



where the angle = corresponds to an electron emission within the reaction plane. 



W^--(p) =c4sin^^ 



(25) 



W """(p) =h + h2 sin^ e cos^ + 64 sin^ 9 cos^ 0, 



(26) 
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Figure 2. Angular distributions of the emitted electrons in the two-photon K-sheW 
ionization of hydrogen- like ions by means of circularly, linearly and unpolarized light. 

Results are presented for both, the electric dipole ( ) and the relativistic { — ) 

approximations and for a two-photon energy which is 40 % above the (one-photon) 
ionization threshold. 



3. Results and discussion 

For the calculation of total two-photon ionization cross sections, the electric dipole 
approximation was recently found sufficient for most of the hydrogen-like ions, and not 
just in the low-Z domain (Koval et al 2003). Even for high-Z ions, for example, the total 
cross sections from the dipole approximation do not differ more than about 20 % from 
those of a full relativistic computation, including the contributions from all the higher 
multipoles. Larger deviations, however, can be expected for the angular distribution of 
the emitted electrons which is known to be sensitive to the retardation in the electron- 
photon interaction. As known, for instance, from the radiative recombination of high-Z 
ions, a significant change in the angle-differential cross sections may arise from the higher 
multipoles and may lead to quite sizeable deviations when compared with the dipole 
approximation (Surzhykov et al 2002). 

In this contribution, therefore, we have analyzed both, the electric dipole and 
the exact relativistic treatment from Eq. (122|) in order to explore the relativistic and 
retardation effects on the angular distributions of the electrons. Detailed computations 
have been carried out, in particular, for the i^-shell ionization of (neutral) hydrogen 
as well as hydrogen-like xenon and uranium. Moreover, to explore the dependence 
of the relativistic effects on the polarization of the incoming light, three cases of the 
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Figure 3. Angular distributions of the electrons emitted in the two-photon if-shell 
ionization of the hydrogen-like uranium U^^+ by means of linear polarized light. 
Distributions are shown for the angles = 0°, 45° and 90° with respect to the reaction 
plane; cf. Figure 1. 

polarization are considered: (i) completely circular polarized, (ii) completely linear 
polarized as well as the case of (iii) unpolarized light. For these three ions and types 
of polarization, Figure 2 displays the angular distributions of the electrons as obtained 
by the dipole approximation ( — ) and the exact relativistic treatment ( — ). While, for 
hydrogen, both approximation virtually yields identical results, the start to differ as the 
nuclear charge Z is increased. Instead of a symmetrical emission with respect to the 
polar angle 6 = 90°, then the emission occurs predominantly into forward direction, an 
effect which is best seen for hydrogen-like U^^~^ ions. We therefore find, that the non- 
dipole terms give first of all rise to an asymmetrical shift in the angular distribution 
of the electrons which could be observed in experiment. The maxima in the (angle- 
differential) cross sections, on the other hand, are less affected and deviate, even for 
hydrogen-like uranium, less than a factor of 2. 

In Figure 2, all angular distributions are shown as function of the polar angle 6, 
i.e. with respect to the incoming photon beam. As discussed above, this dependence 
of the differential cross sections, ^ = ^(^) , can be the only one for circular and 
unpolarized light for which the electron emission must be axially symmetric. For linear 
polarized light, in contrast, the emission of the electrons will depend on both, the polar 
angle 6 and the azimuthal angle 0. For this polarization. Figure 2 only displays the 
angular distributions within the reaction plane, i.e. at = 0°. To explore, in addition, 
also the 0-dependence of the two-photon ionization by linear polarized light explicitly. 
Figure 3 shows the corresponding angular distributions ^{0, 0) for the three particular 
angles = 0°, 45° and 90° with respect to the reaction plane; here, the left inlet (0 
= 0°) is the same as shown in Figure 2 in the middle column for U^^^ ions. Again, 
the results from the electric dipole approximation are compared with those from a fully 
relativistic computation. As seen from Figure 3, the most pronounced effect of the 
higher multipoles arise for an electron emission in a plane, which is perpendicular to the 
photon polarization vector (0 = 90°). In such a — perpendicular — geometry of the 
experiment, the cross sections from the exact treatment show strong forward emission 
of the photoelectrons while the electric dipole approximation (12^ . in contrast, results 
in a completely isotropic emission, if seen as function of the polar angle 6. 
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Figure 4. Angular distributions of the electrons emitted in the two-photon ii'-shell 

ionization of the hydrogcn-likc uranium U^^^ by circular polarized light with different 
degrees of polarizations Pc = 0, 0.5, 0.7. and 1. 

Until now, we considered the two-photon ionization of hydrogen-hke ions by either 
completely polarized (linear: Pl = 1; circular Pc = 1) or unpolarized light {Pl = Pc = 
0). In most experimental investigations on two- (and multi-) photon processes, however, 
the incident radiation is typically polarized with some given degree of polarization 
< Pc, Pl < 1- Apart from the type of the polarization of the incoming light, 
therefore, we shall study also how the angular distributions depend on the degree or 
polarization. Figure 4, for instance, displays the angular distribution from the shell 
ionization of hydrogen-like U^^"'" ions by means of circular polarized light with a degree 
of polarization Pc = 0.0 (unpolarized case), 0.5, 0.9 and 1.0. As seen from this 
figure, the probability for an electron emission increases at angles around 6 = 60° but 
decreases (towards zero) in forward and backward direction as the degree of polarization 
is increased. In particular the behaviour near ^ = 0° and 6 = 180° can be easily explained 
if we consider the conservation of momentum in the overall system. Since, for completely 
circular polarized light the (total) angular momentum of photons on the quantization 
axis (which is chosen along the photon momenta k) becomes Ai -|- A2 = ± 2, it obviously 
can not be compensated — in the final state — if the electron is emitted parallel (or 
antiparallcl) to the incoming light, and hence its spin projection is /i/ = = ±1/2. 
For unpolarized light, in contrast, the photons may have different helicities and, hence, 
their angular momentum Ai + A2 ± will be conserved under a forward and backward 
non spin-flip electron emission. 

4. Summary 

In this paper, the two-photon ionization of hydrogen-like ions has been studied in the 
framework of second-order perturbation theory and the rclativistic description of the 
electron and photon fields. That is, exact Dirac bound and continuum wave functions 
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were applied for the description of the electron to reveal the importance of relativity 
on the angular distributions of the emitted electrons. Moreover, relativistic Coulomb- 
Green's functions are used to perform the summation over the complete Dirac spectrum 
as needed in second-order perturbation theory. 

To understand the angular distributions of the emitted photoelectron and, in 
particular, the influence of the polarization of the light on this emission, density matrix 
theory has been utilized to 'combine' the two-photon transition amplitudes in a proper 
way. Calculations are carried out for the K-sheW ionization of the three (hydrogen-like) 
ions H, Xe^^+ and U^^+. From the angular distribution of the electrons for different 
types (linear, circular, unpolarized) and degrees of polarization (i.e. in going from the 
completely polarized to unpolarized light), it is clearly seen that the angular emission 
depends much more sensitive on the contributions from higher multipoles than the total 
cross sections. Two rather pronounced effects, for example, concern the (asymmetrical) 
forward emission of the electrons as well as a significant change in the electron emission 
for linear polarized light, if the electrons are observed perpendicular to the reaction plane 
[cf. Figure 4]. Both effects are enhanced if the nuclear charge of the ions is increased. 

An even stronger influence from the non-dipole terms (of the radiation field) is 
expected for the spin-polarization of the photoelectrons. Similar as in the present 
investigation, density matrix theory provides a very suitable tool for such polarization 
studies. A detailed analysis of the polarization of the photoelectrons, emitted in the 
two-photon ionization of hydrogen-like ions, is currently under work. 



Appendix A. Photon spin density matrix 

A pure (i.e. completely polarized) state of the photon can be characterized in terms 
of a polarization unit vector u which always points perpendicular to the (asymptotic) 
photon momentum k. Of course, this polarization vector, u, can be rewritten by means 
of any two (linear independent) basis vectors such as the circular polarization vectors 
u-ti which are (also) perpendicular to the wave vector k and which, for u+i respective 
u_i, are associated with right- and left-circular polarized photons (BlumSl). In such a 
basis, the unit vector for the linear polarization of the light can be written as 

u(x) = ^ ( e-^'^ u+i + e^'^ u_i) , (A.l) 

where x is the angle between u(x) and the x-z plane. 

While a description of the polarization of the light in terms of either the circular 
polarization vectors uj-i or the linear polarization vector flA.l|) is appropriate for 



completely polarized light, it is not sufficient to deal with an ensemble of photons which 
have different polarization. Such a — mixed — state of the light is then better described 
in terms of the spin-density matrix. Since the photon (with spin S = 1) has only two 
allowed spin (or helicity) states |kA) , A = ±1, the spin-density matrix of the photon 
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is a 2 X 2 matrix and, hence, can be parameterized by three (real) parameters: 

I f 1 + Pc Pl e-'^^ \ 
kA|p,|kA) = - , (A.2) 

/ ^ \ Pl e^^^ I - Pc ) 

where < Pl < 1 and — 1 < Pc < 1 denote the degree of hnear and circular polarization, 
respectively. The angle Xi nioreover, represents the direction of the maximal linear 
polarization of the light. 

Of course, the choice of the parameters Pi, Pc and x is not unique and many other 
— equivalent — sets of three real parameters could be applied to characterize the photon 
spin density matrix ( ]A.2I) . In the analysis of experimental data, for instance, one often 
uses the three Stokes parameters to describe the polarization of radiation. The Stokes 
parameters can easily be expressed in terms of the (two) degrees of polarization, Pl and 
Pc, and the angle x as: 

Pl = Pl cos 2x, P2 = Pl sin 2x, P3 = Pc ■ (A.3) 

The use of the Stokes parameters leads to the familiar form of the spin density matrix 
(Balashov et al 2000) 



1 



1 + P3 Pl - tP2 



kA |/3,|kA ) = - I _ _ _ I . (A.4) 



Pl + zP2 1 - P3 
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